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Amyotrophic lateral sclerosis (ALS) is a devastating, rapidly progressing disease of the central nervous 
system that is characterized by motor neuron degeneration in the brain stem and the spinal cord. We 
employed time of flight secondary ion mass spectrometry (ToF-SIMS) to profile spatial lipid- and 
metabolite- regulations in post mortem human spinal cord tissue from ALS patients to investigate chemical 
markers of ALS pathogenesis. ToF-SIMS scans and multivariate analysis of image and spectral data were 
performed on thoracic human spinal cord sections. Multivariate statistics of the image data allowed 
delineation of anatomical regions of interest based on their chemical identity. Spectral data extracted from 
these regions were compared using two different approaches for multivariate statistics, for investigating 
ALS related lipid and metabolite changes. The results show a significant decrease for cholesterol, 
triglycerides, and vitamin E in the ventral horn of ALS samples, which is presumably a consequence of motor 
neuron degeneration. Conversely, the biogenic mediator lipid lysophosphatidylcholine and its fragments 
were increased in ALS ventral spinal cord, pointing towards neuroinflammatory mechanisms associated 
with neuronal cell death. ToF-SIMS imaging is a promising approach for chemical histology and pathology 
for investigating the subcellular mechanisms underlying motor neuron degeneration in amyotrophic lateral 
sclerosis. 



Amyotrophic Lateral Sclerosis (ALS) is a very progressive, fatal disease that is characterized by irreversible 
degeneration of motor neurons in the spinal cord, brain stem, and cortex resulting in increasing muscle 
weakness and muscle atrophy^ ALS has a prevalence of 2-4 cases per 100 000 a year with a common onset 
age between age 40 and 60 and an average life expectancy three to five years after diagnosis^. Most ALS cases are 
idiopathic in nature; however, 10% have been found to be caused by genetic factors^. The main mechanisms 
underlying ALS are still unknown, although a couple of findings concerning molecular ALS pathology have been 
reported mainly focusing on genetic causes^"^. The limited knowledge of the disease results in further drawbacks 
such as the inaccuracy and time delay for diagnosis, which highlights the need for further investigation of the 
disease on a molecular level. 

The limited understanding of molecular ALS pathobiology relates directly to the lack of appropriate analytical 
technologies that feature the necessary sensitivity and specificity. Studying complex biological and pathological 
processes requires acquisition of molecular images, since these processes involve the translocation and interaction 
of a wide range of chemical species. 

Mass spectrometry based molecular imaging (IMS) is an emerging technique for probing the molecular 
architecture of biological tissue in an anatomical context^. In contrast to common molecular biology techniques, 
IMS does not require any a priori knowledge of the potential target species and allows direct and comprehensive 
identification and spatial profiling of biomolecular species in situ, while maintaining high molecular specifi- 
city^"^°. A major advantage of ToF-SIMS based imaging MS is the possibility to robustly acquire high- resolution 
images at the cellular length scales (1-5 jim), without compromising sample treatment procedures^^"^^ These 
features make SIMS a powerful technology for spatially profiling lipids and metabolites at the single cell leveP^. 
This is of particular relevance for studying complex and heterogeneous samples, such as neuronal tissue and cells, 
which arguably constitute the most complex and least understood systems in the human body^'^^. ToF-SIMS 
based molecular imaging represents therefore a promising approach for investigating spinal cord biochemistry. 
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Figure 1 | Image Data Analysis for Tissue Segmentation. Multivariate statistical analysis of imaging MS data by means of maximum autocorrelation 
factor analysis (MAF) reveals anatomical regions of human spinal cord including grey and white matter in both control and ALS samples. Composition 
images of corresponding factors (i.e. components, MA- factors) allow identification of anatomical regions of interest. Each pixel intensity represents the 
respective score (eigenvalue, see scale). Scale bar = 4 mm. 



Several studies on spinal cord using various imaging mass spec- 
trometry techniques have been reported that have provided import- 
ant data on spinal cord pathology^^"^^. In a recent study, we probed 
the native distribution of lipids and metabolites in human spinal cord 
using ToF-SIMS for high spatial resolution imaging^^. In addition, 
different strategies for multivariate analysis of whole image datasets 
for tissue segmentation were evaluated that in turn facilitate iden- 
tification of anatomical regions of interest based on their chemical 
identity. 

Here we report the use of ToF-IMS based analysis of post mortem 
human spinal cord tissue samples in conjunction with multivariate 
data analysis (MVA) to identify ALS related changes in spinal cord 
lipid chemistry and metabolism in situ. This work should gain fur- 
ther insight in molecular mechanisms underlying ALS, particularly 
with respect to lower motor neuron degeneration. 

Results 

Multivariate image analysis for identification of anatomical re- 
gions of interest in human spinal cord. For each individual (n = 5; 2 
control/3 ALS), subsets of two consecutive tissue sections were 
analyzed with ToF-SIMS in both positive and negative ion mode. 
The datasets were calibrated internally and reconstructed with a 
common mass interval list (binning) for data acquired in both 
positive- and negative ion mode. This resulted in a total of 20 
imaging datasets (10 pos.- and 10 neg. mode), each at least 3 GB 
in size. Imaging data from the 20 datasets were subjected to 
unsupervised multivariate statistical image analysis by means of 
maximum autocorrelation factor analysis (MAF)^^'^^. This resulted 



in generation of eigenvalue (score) images that depict the chemical 
variance captured in the respective factors of the MAF analysis 
(Figure 1). This image analysis approach revealed characteristic 
outlines of the grey and white matter (Figure 1). The molecular 
species that outline the histological regions and cellular environ- 
ments can then be deduced from the loadings (eigenvectors) and 
provide further information on local lipid biochemistry and 
metabolism. 

Multivariate analysis of region of interest data allows discrimina- 
tion of ALS and control groups. Two different multivariate 
approaches, partial least square discriminant analysis (PLS-DA) 
and "Statistical Analysis of Microarray Data" (SAM), were 
employed for statistical analysis of spectral data extracted from 
anatomical regions of interest (ROI). These histological regions 
including the whole grey matter (GM) as well as the ventral horn 
of the grey matter (VH) were assigned in accordance with the image 
analysis data (Figure 2A). The ROI data of the GM and the VH of 
each patient were normalized to the total ion count (TIC), log2 
converted, and different ionization modes were pooled for each 
patient. This gave two sample/variable matrices for both the GM 
ROI and the VH ROI that were then evaluated by means of 
multivariate analysis using two different strategies. Here, as a first 
approach for MVA, multivariate analysis by means of PLS-DA was 
performed. PLS-DA is a linear regression extension of principal 
component analysis (PGA) and well suited for comparing two 
different sample groups. PLS-DA was performed after mean center- 
ing and Poisson scaling and resulted in distinct discrimination of the 
two groups (Gontrol and ALS) both for the grey matter region and 





Figure 2 | Multivariate statistics of ROI spectral data using PLS DA. MS spectral data from previously outlined anatomical regions (A), including the 
grey matter (GM) (B) and ventral horn (VH) (C), were subjected to multivariate analysis using PLS DA. (B,C) PLS DA score plots demonstrate grouping 
of ALS samples (green) and controls (red). 
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Figure 3 | Statistical Analysis of Spectral Data using SAM. Spectral data of previously outlined regions as revealed by image analysis were subjected to 
statistical analysis of microarray data (SAM). Spectra acquired in positive and negative mode for the GM- and VH-areas (A) were analyzed. (B, C) SAM 
plots illustrate relative changes of ALS versus control in the respective areas as revealed by unpaired t-statistics. The x-axis illustrates the degree of 
significance, while the y-axis refers to the observed degree of change. Red illustrates a relative increase in ALS compared to control, while green illustrates a 
relative decrease in ALS to control. The borders indicate the significance threshold (p < 0.05). 



even more pronounced for the ventral horn data (Figure 2B, C). 
From the corresponding loading plots the individual variables {m/z 
values) were deduced that have the greatest impact on the score and 
thereby separation (SI Figure lA, B). 

As a second approach, the "Statistical Analysis of Microarray data" 
tooP^ was used for grouped unpaired t-statistics of the ALS and con- 
trol group data for each ROI (GM and VH) (Figure 3A). For the whole 
grey matter region the results show significant intensity changes (p < 
0.05) for in total 21 m/z values. Here, four mass peaks were signifi- 
cantly increased in ALS (red), while a number of 17 species were 
higher in the control group (green. Figure 3B). For the ventral horn 
region; however, the data show significant changes (p < 0.05) for 
175 m/z values that were found either increased (130) or decreased 
(45) in the ALS samples compared to the controls (Figure 3G). 

For verification of these fmdings, SAM results for the respective 
ROI were compared to the results obtained from PLS-DA. Here, the 
PLS-DA derived variables (loadings) that exhibited the greatest con- 
tribution to the corresponding score of the individual samples were 
investigated. Only, m/z compounds that were found significant in 
SAM and displayed a large PLS-DA loading value were considered. 
For the grey matter region all 21 compounds that were found sig- 
nificantly changed in between ALS and controls were also found to 
have a large loading value in PLS-DA. In contrast, of the 175 signifi- 
cant observed in SAM analysis for the VH region, only 40 were found 
to have a large loading value. These were then further verified using 
non-parametric unpaired statistics (Kruskal-Wallis). This resulted in 
a number of 38 compounds that were found significantly changes in 
the ventral horn in between ALS and control samples. A number of 
15 peaks were increased and 23 peaks were decreased in ALS ventral 
horn as compared to controls. Here, for the decreased compounds all 
peaks were detected in negative mode, while for increased com- 
pounds 7 were detected in positive and 8 in negative mode. 

Regional changes of spinal cord lipids associated with ALS. The 

most prominent signals were tentatively assigned based on accurate 
mass matching and comparison to values reported in the litera- 
ture^^"^^ as well as by matching the experimental data to predicted 
isotopic patterns. Multivariate analysis of the ROI spectral data 
resulted in 6 molecular species that were found decreased in the 
grey matter as well as in the ventral horn of ALS spinal cord. 
Five of these peaks could be assigned to triglyceride species (TG, 
Figure 4A)^^'^^, while one compound {m/z 446.3) could not be 



identified. Further mass peak results obtained for the whole GM 
ROI could not be verified by manual inspection of the single ion 
images. By confining comparative analysis to the VH region more 
species were found significantly changed in between ALS and 
controls. Among the most prominent compounds that displayed an 
ALS associated decrease, cholesterol (C27H45O: [M — H]" 385.3) and 
vitamin E (C29H49O2: [M - H]" 429.3) could be identified"". For 
both species a significant decrease was observed in the ventral horn 
region of controls as compared to ALS samples (Figure 4B). 

Conversely, a number of compounds were found increased in the 
ventral horn of ALS tissue as compared to controls. This included 
phosphatidylcholine (PC) -lipid fragments (TME, trimethylethyli- 
mine: m/z 86.1, [C5H12N+]; PC-headgroup: m/z 184.1 [C5H15 
PO4N+]) as well as intact lysophosphatidylcholine (LPC (16:0), 
m/z 496.5 [M + H]^) (Figure 5A). Furthermore, protein associated 
signals (m/z 26.01, [CN"]; m/z 80.96 [NaKCN+]) and hydrogensul- 
fite [HS03~] were found elevated in the ventral horn of ALS spinal 
cord samples (Figure 5B). 

Discussion 

The spinal cord constitutes the interface of the peripheral and the 
central nervous system comprising essential circuits in motor control 
as well as somatosensory function. Anatomically, the tissue is gen- 
erally separated into the white matter and grey matter. Nervous 
tissue is characterized by its high complexity and diversity. This will 
eventually be reflected in the spatial intensity distribution pattern of 
biochemical species and in turn affect signal intensity values based on 
changes in the biochemical background. Therefore comparison 
between heterogeneous tissue specimens requires identification 
and annotation of distinct anatomical regions so that histological- 
and biochemical similar areas can be compared in between different 
individuals. For comparative analysis, we employed a ToF-SIMS 
imaging data guided approach for tissue segmentation of spinal cord 
samples. The aim was to outline anatomical regions of interest based 
on their chemical identity. This approach was based on multivariate 
analysis of whole scan image data using maximum autocorrelation 
factor analysis. The suitabiHty of this approach was evaluated in a 
recent study by our group, where we compared the performance of 
PGA and MAF for image analysis with respect to contrast and their 
ability to identify anatomically relevant regions of interest^^. Here, a 
MAF based image analysis guided strategy for ROI assignment was 
found to be a robust and straightforward approach for tissue seg- 
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Figure 4 | Single Ion Images of Chemical Species decreased in ALS. Multivariate statistics of ROI spectral data using SAM and PLS DA resulted in 
detection of multiple molecular species displaying ALS associated decrease in the grey matter. Only compounds observed in negative mode displayed a 
decrease in ALS spinal cord. (A) Here, a decrease throughout the whole grey matter (arrows) was observed for various triglycerides (TG (46 : 4), [M — H]~ 
769.6; TG (46 : 3), [M - H] " 771.6; TG (50 : 1), [M - H] " 831.7). (B) Similarly, an unidentified peak at m/z 446.3 was found decreased in both ventral 
horn and whole grey matter (arrows). In contrast, cholesterol (C27H45O, [M — H]~ 385.3) and vitamin E (C29H49O2, [M — H]^ 429.3) were found 
consistently decreased in the ventral horn of ALS spinal cord but not in the whole grey matter. Scale bar = 4 mm. 



mentation^^. MAF employs an implemented correction to the noise 
covariance matrix as achieved by an iterative comparison to adjacent 
pixel spectra. In contrast to PCA, this approach is therefore not 
dependent on data preprocessing (e.g. scaling, mean centering) 
and thus provides better contrast and sensitivity that allow catching 
slight intra tissue heterogeneities, although the technique is compu- 
tational demanding. This highlights an important consideration for 
bulk- image data analysis with respect to sample size and data hand- 
ling. Given the large data file size obtained from the SIMS and other 
IMS experiments (up to 10 GB), tradeoffs with respect to pixel reso- 
lution and number of peak assignments (bins) have to be made. 
Given the fact that most commercial solutions for image data ana- 
lysis, as provided by the manufacturers, are still limited to 32 -bit 
systems, the complexity of downstream data analysis approaches 
have to be considered in advance. In the present study, MAF facili- 
tated exact annotation of these regions and was employed for sub- 
sequent extraction of ROI spectral data for comparative analysis of 
the localized chemistry in ALS tissue compared to controls. Using 
PLS-DA as well as SAM based statistics, the most significant changes 



were observed by confined comparative analysis of the ventral horn 
region (Figure 2 and 3). These observations further demonstrate that 
targeted ROI data analysis is mandatory to reveal local, disease 
related, biochemical changes. Indeed, ALS is characterized by degen- 
eration of spinal cord motor neurons in the ventral horn. Spatially 
confined changes in lipid and metabolite expression are easily con- 
voluted, if the ROI is assigned too widely, which further illustrates the 
relevance for targeted imaging data analysis when analyzing hetero- 
geneous biological tissue. 

A decrease of triglyceride species was observed throughout the 
whole grey matter (Figure 4A). This is consistent with previously 
reported results showing an ALS associated increase in clearance of 
peripheral lipids, including triglycerides^^. In that work, an altered 
energy metabolism was observed in experimental ALS (mSODl mice) 
that resulted in increased TG metabolism. This shift in energy meta- 
bolism relates to an increased demand of nutrients in muscle tissue as 
a consequence of dysfunctional skeletal muscle metabolism^°. 

While 21 m/z peaks were found statistically changed between ALS 
and control spinal cords for the whole GM region, they could not be 
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Figure 5 | Single Ion Images of Chemical Species elevated in ALS spinal cord. Multivariate statistics of ROI spectral data using SAM and PLS DA 
resulted in detection of several molecular species displaying ALS associated increase in the grey matter. (A) These included lysophosphatidylcholine (LPC 
( 16 : 0), [M + H] ^ 496.5, arrow) as well as its fragments, PC-headgroup (PC) and trimethylethylimine (TME), observed in pos. ion mode. (B) In addition, 
protein fragments (m/z 26.01, [CN~] and m/z 87.9, [NaKCN^] ) as well as hydrogensulfite (m/z 80.96, [HS03~] ) were found elevated in the ventral horn 
of ALS spinal cord (arrow). Scale bar = 4 mm. 



verified by manual inspection of the single ion images. This in turn 
highlights the importance of strict validation of MVA results by 
subsequent inspection of the imaging raw data. As mentioned prev- 
iously, further refinement of the GM-ROI to the ventral horn yielded 
a significantly higher number of ALS associated molecular species 
(Figure 3C). Here, the data show that cholesterol (C27H45O: [M — 
H]- 385.3) and vitamin E (C29H49O2: [M - H]" 429.3) were 
increased in the ventral horn region of controls as compared to 
ALS samples (Figure 4B). Cholesterol and vitamin E are major con- 
stituents of the myelin sheath surrounding ascending and descend- 
ing axons that constitute the white matter but also innervate the grey 
matter regions. These results are consistent with the fact that ALS 
is characterized by neuronal degeneration in the ventral horn. 
Decreased levels of the myelin sheath constituents correspond to 
decreased innervation as a result of neuronal cell death. Moreover, 
the degree of decreased white matter compounds may serve as a 
neuropathological hallmark of neurodegenerative mechanisms 
underlying ALS. 

An ALS associated increase of lipids and metabolites was also 
observed and verified in the ventral horn (Figure 5). Here, 
SAM and PLS-DA revealed an increase in PC-lipid fragments 



(PC-headgroup and TME, Figure 5 A) in the ventral horn of the grey 
matter in ALS patients (Figure 5A). Similarly, LPC (16:0) was sig- 
nificantly increased in the same region in ALS spinal cord tissue, 
while no other phospholipids displayed any significant change 
(Figure 5A). PC-headgroup and TME are most likely fragments of 
LPC resulting from in- source fragmentation. The signal can however 
also originate from sphingomyelin^^ A mass peak at m/z 224.1 cor- 
responding to the PC-lipid specific fragment C8Hi9PN04^ (PC- 
headgroup + C3H4) displayed similar distribution pattern and regu- 
lation in intensity as the other PC fragments (SI Figure 2)^\ This 
verifies that the PC attributed mass peaks {m/z 184.1; m/z 86.1) do 
indeed originate from PC-lipids, rather than sphingomyelin. 

LPC is a precursor for lysophosphatitic acid (LP A), a biogenic 
lipid. This lipid mediator is synthesized from LPC by autotaxin 
(ATX) that was demonstrated to possess lysophospholipase D activ- 
ity^^. Interestingly, in cerebrospinal fluid (CSF) the levels of LPC and 
LP A are very low despite significant amounts of ATX, which indi- 
cates that LPC production would be a rate-limiting step for LPA 
production^l LPA production and activity in the spinal cord has 
been described before, where LPA was found to exert its activity 
through LPAl -receptor signaling, which initiates neuropathic pain 
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and underlying mechanisms including demyelination^^. Elevated 
LPA levels are a consequence of serum leakage through tissue 
damage but also result from abnormally elevated LPA biosynthesis^^ 
LPA/LPAl -receptor signaling in the spinal cord was found to initiate 
e.g. neuropathic pain through different demyelination pathways^l 
These include phosphoinositol-bi/tri-phosphate (PIP2/PIP3) second 
messenger mediated genomic regulation of myelin protein express- 
ion^^ as well as inositoltriphosphate (IP3) mediated myelin degrada- 
tion by calpaein^^. This is consistent with the observations presented 
here, where decreased cholesterol and vitamin E values indicate 
demyelination. In addition to these observations, a significant 
increase of protein specific fragments {m/z 26.01, [CN"]) was 
observed in the ventral of the ALS spinal cord (Figure 5B). Similar 
to other neurodegenerative diseases, ALS is characterized by pro- 
gressive accumulation and aggregation of misfolded proteins into 
insoluble protein aggregates. Increased protein signals therefore 
indicate impaired protein turnover. Again, this is consistent with 
previously reported observations, where a diminished proteasome 
activity has been observed in the ventral horn of ALS patients^^'^^. We 
also observed increased levels of hydrogen sulfite {m/z 80.9, 
[HSOs"]) in the ventral horn of ALS spinal cord. Sulfite is a meta- 
bolite of cysteine, which accumulates when one or more of the 
enzymes needed for glutathione (GSH) synthesis are defective^°. 
Sulfite exerts its neurotoxic potential through direct inhibition of 
glutamate dehydrogenase thereby impacting ATP biosynthesis^^. 
The present observations on spinal cord sulfite accumulation indi- 
cate cysteine and sulfite mediated neurotoxicity as a consequence of 
disrupted GSH metabolism, as has been discussed earlier^\ 

Conclusions 

The potential of ToF-SIMS imaging for interrogating the biochem- 
ical architecture of the human spinal cord in ALS is demonstrated. 
Multivariate statistics used for image data guided tissue segmenta- 
tion facilitate accurate outlining of histological regions of interest 
based on their chemical identity. Spectral data analysis of these 
ROI's revealed spatially confined changes in lipid and metabolite 
regulation associated with ALS pathogenesis. ALS implicated 
decrease of cholesterol and vitamin E in the ventral horn was 
observed presumably associated with demyelination and motor neu- 
ron degeneration. In addition, ventral horn triglyceride clearance 
indicates ALS related changes in energy metabolism, while an accu- 
mulation of hydrogen sulfite and LPC in the ventral horn point to 
metabolic neurotoxicity and biogenic lipid induced neuroinflamma- 
tion, respectively. The data presented reveal novel targets for investi- 
gating ALS pathology in spinal cord that have not been presented 
before. ToF-SIMS can be used to reveal disturbances in the neuro- 
chemical interplay affecting neuronal functionality and integrity and 
suggests a novel approach for future investigations concerning motor 
neuron degeneration. 

Methods 

Chemicals and Reagents. All chemicals were of pro-analysis grade and purchased 
from Sigma Aldrich (St. Louis, USA). TissueTek optimial cutting tool (OCT) was 
purchased from Sakura Finetek (AJ Alphen aan den Rijn,The Netherlands). Water 
was purified with a Milli-Q (Millipore, Bedford, MA, USA) purification system. 

Clinical Samples. Post-mortem spinal cord tissue samples of thoracic level from two 
controls and three ALS patients were used in the study. The cause of death in the two 
control subjects was cardiac failure in one case and myocardial infarction in the other. 
At autopsy, the spinal cords were dissected, removed from the dura mater, and cut 
into 5-mm sections. The tissue was then immediately frozen on a metal plate 
maintained in liquid nitrogen and stored at — 72°C. Autopsy was performed with 
written approval from the family in accordance with Swedish rules and regulations. 

Sample Preparation. Thoracic spinal cord sections (12 |im) were collected using a 
cryostate microtome (Thermo Scientific, Waltham, MA) and thaw- mounted on 
conductive glass slides (indium tin oxide, ITO, Bruker Daltonics, Bremen, Germany). 
The thaw mounted sections were dried under vacuum for 10 min and stored at 
— 80°C until further use. For SIMS experiments, 2 consecutive sections from each 
subject were analyzed. 



ToF-SIMS Analysis. All ToF SIMS experiments were performed on an lONTOF V 
instrument (lONTOF GmbH, Miinster, Germany) equipped with a Bis"^ cluster ion 
gun as primary ion source. For whole section analysis, data were acquired in high 
current bunched mode (HBC)^^ with a primary ion current of 0.25 pA at 25 keV and 
a maximum ion dose was 1.5 10^/cm^. Scans were acquired in positive and negative 
ion mode, using the "stage scan macro raster function" with 10 shots per pixel on 0.4 
X 0.4 mm areas with 200 measurements per mm resulting in a pixel resolution of 
5 |im. This acquisition mode comprises stepwise acquisition of 0.4 X 0.4 mm 
patches. Each patch was acquired in electrostatic raster mode where the beam was 
moved over the field of view. The whole tissue area was analyzed by stepwise 
acquisition of these patches as achieved by stepwise movement of the sample stage. 
The mass resolution in high current bunch mode was about M/AM = 5000. 

Data Processing. All data were acquired and processed with the Surface Lab software 
(v. 6.3 lONTOF). All spectra were calibrated internally to signals of [C]"^, [CH]"^, 
[CHa]^, [CH3]+, [C5Hi5PN04]+ and [C27H45]+ that were used as calibration points 
in positive mode and [C] [CH] [O] ", [C16H31O2] ", [C18H35O2] " and [C27H45O] " 
in negative ion mode. Two mass interval lists were created (one for each ion polarity), 
which contained the list of m/z values to be included in multivariate analysis. Mass 
interval lists were created by peak search of all individual samples according to the 
following search parameters: S/N > 3, width 0.8 Da. The assigned mass peaks were 
collected in the same mass interval list and redundant peak assignments removed. 
Image data were exported into the *.bif6 format and spectral data were exported 
directly from the peak statistics tab into Excel (v. 2010). 

Data Analysis. For statistical analysis two approaches were employed comprising 
intact image data analysis as well as statistical analysis of spectral data of distinct 
regions of interests. For image analysis, *.bif6 files were loaded into PLS-Toolbox (v. 
7.02, Eigenvector Research Inc., Wenatchee, WA) running in MatLab (v. 12a, The 
Math Works Inc., Natick, MA). Multivariate image analysis was performed by means 
of maximum autocorrelation factor analysis (MAF). For comparative statistics on 
spectral data, relevant anatomical regions of interest including the grey matter as well 
as the ventral horn of the grey matter were outlined followed by reconstruction of the 
respective ROI using the corresponding feature in the Surface Lab software 
(lONTOF). Two different approaches were employed for comparative, statistical 
analysis of ROI spectral data. These included 1) partial least square discriminant 
analysis (PLS-DA) as well as 2) "Statistical Analysis of Microarray data" (SAM)^*^. 
Peak area values were log2 normalized followed by SAM based grouped unpaired 
statistical analysis (t- statistics) using the corresponding macro in ExceP*^. PLS-DA 
was performed on the spectral data using PLS-Toolbox (Eigenvector Research). Here, 
data were mean centered and Poisson scaled. Further analysis of single mass peaks 
and comparisons between groups was performed with non-parametric ANOVA 
(Kuskal- Wallis) with 95% confidence interval. Identification of mass peaks was done 
using tentative assignment by comparing the accurate mass to previously reported 
values in the literature^^ as well as matching of isotopic distribution. 
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